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' We study dynamical moduli stabilization driven by gaugino condensation in supergravity. In 

0\ . the presence of background radiation, there exists a region of initial conditions leading to 

' successful stabilization. We point out that most of the allowed region corresponds to initial 

Hubble rate H close to the scale of condensation A. which is the natural cutoff of the effective 
' theory. We first show that including the condensate dynamics sets a strong bound on the 

, initial conditions. We then find that (complete) decoupling of the condensate happens at H 

' about two orders of magnitude below A. This bound implies that in the usual scenario with 

the condensate integrated out, only the vicinity of the minimum leads to stabilization. Finally, 
we discuss the effects of thermal corrections. 



1 Introduction 



Higher dimensional supersymmetric theories such as string theory are good candidates to explain 
the origin of four-dimensional physics at low energies. In general, compactification of the addi- 
tional dimensions yields a large vacuum degeneracy, parametrized in terms of the so-called moduli 
fields. In order to reproduce our observable universe and its particle content, consistent higher di- 
mensional theories must offer some mechanisms to fix the mass and the vacuum expectation value 
(vev) of the moduli. In type IIB string compactifications, nontrivial background fluxes stabilize 
most of the geometrical moduli [T], but leave some other flat directions unlifted. In particular, the 
overall volume parametrized by a Kahler modulus a is not fixed and hence requires the inclusion 
of nonperturbative effects [2H4] . 

Gaugino condensation [36] is perhaps the best understood and most successful such mechanism 
capable of providing masses to moduli. It is a key ingredient in racetrack models [7], models with 
Kahler stabilization [8], and the Kachru-Kallosh-Linde-Trivedi (KKLT) scenario [9]. The success 
of these models relies on the assumption that the condensate forms at a high scale, inducing an 
effective scalar potential for moduli. However, it is known that dynamical stabilization of moduli is 
fragile and depends on the particulars of the universe evolution |10| . If the cosmological expansion 
is dominated by some energy source other than the moduli, such as radiation or vacuum energy, 
the initial conditions leading to stabilization may be less restricted than in a modulus dominated 
universe. This statement has been confirmed in matter and radiation dominated scenarios |11H13) 
including finite temperature corrections |14| . 

Despite this progress, there is no reason to believe that only moduli dynamics will be affected 
by the evolution of the universe. Indeed, it has been recently shown that in a universe undergoing 
fast expansion, the condensate itself may be destabilized leading to unsuccessful moduli stabiliza- 
tion |15| . Thus, also the condensate is altered by the presence of additional sources of energy. 
Whether gaugino condensation occurs depends on the dynamics of the gaugino pairs, which can be 
described in terms of a complex (super)field u in the Veneziano-Yankielowicz (VY) approach [5]. 
In the present work we take this view seriously and study the consequences of the interaction 
between the condensate and background radiation upon moduli stabilization. The main questions 
we will be concerned with are: 

• when can the condensate field u be safely integrated out without introducing inconsistencies 
in the theory? 

• for which initial conditions does gaugino condensation lead dynamically to moduli stabiliza- 
tion? 

We shall show that for most configurations that lead to moduli stabilization, the condensate cannot 
be integrated out and its dynamics must be considered. Although, for illustrative purposes, we 
shall focus on the KKLT model, it is straightforward to extend our findings to any model of moduli 
stabilization driven by gaugino condensation, where similar results are expected. 

2 Single field dynamics 

As mentioned above, we consider the KKLT model [9] for moduli stabilization in the context of 
type IIB string theory. After fixing the complex structure moduli and the dilaton, the effective low 
energy Af = 1 four-dimensional supergravity theory only contains the volume modulus a, which is 
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stabilized by gaugino condensation0 More specifically, an SU(A^) gauge theory localized on D7- 
branes wrapped around the compactification volume has a field-dependent coupling = 1/cJr . 
If this sector is asymptotically free, it undergoes gaugino condensation below a dynamical scale 
A. Using the VY approach [5], it has been shown |17| that the only Kahler and superpotential 
consistent with all symmetries are 

K = -31nf a + a- ^) , (2.1) 



3 

W = {a + 2c\iiu) , 

where the chiral superfield oc Tr W°'W°' admits the condensate A"A° as lowest component. The 
constant c = 3A^/ (iGvr^) is the one-loop beta function coefficient. 
In supergravity, the scalar potential is given by 



K'WiWD.W - 3|T^|^ 



(2.2) 



where K'^ = (K-.j)-' is the inverse Kahler metric, and DiW = Wi + KiW is the usual covariant 
derivative defined over the Kahler manifold. In the limit of small li <C 1, the dominant contribution 
to the potential is 

Vf ~ e^|D„VFpif"", (2.3) 
which has as extrema = 0, — and Wuu — 0, or 

n = , u = Umin = exp {^Yc ~ ^) ' " ^ ^ ~ ^) ■ 

The first solution W = u = Q formally corresponds to a supersymmetric chirally invariant vacuum. 
However, the existence of such a vacuum is inconsistent |18[I19|. Arguably, this minimum can 
be interpreted as an unstable nonsupersymmetric state |20], where the VY potential cannot be 
trusted. The other two solutions correspond to a minimum and a barrier separating it from the 
unphysical state at the origin. Note that they are close in field space. The minimum timin is the 
VY-solution corresponding to gaugino condensation. Thus, we define the condensation scale 

A = Umin • (2.5) 

The mass of the field u around the minimum is of order A, and therefore the field can be integrated 
out at energies below this scale. Plugging the result into (12. ip . expanding in small uu and adding 
a constant Wq arising from integrating out the fiux-stabilized moduli, we recover the familiar 
KKLT potential 

K = -31n(CT + CT) , 

W = Wq + Wnp = Wo + Ae-""" , (2.6) 

and identify A = —1/ (ae), a = 3/ (2c). Note that for reasonable gauge groups, the constant a is 
bigger than one. Also, the inclusion of Wq in (12. ip does not significantly displace the position of 
the extrema (|2^ 



^Note that a similar scenario arises in the weakly coupled heterotic string [16| . where o has to be replaced by 
the dilaton S. 

^From now on, the subscripts r and i denote respectively real and imaginary parts of the fields. 
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Having clarified how gaugino condensation generates an effective potential for the modulus, 
we turn our attention to the single field KKLT model. As is well known, the scalar potential (I2.2p 
for a exhibits a supersymmetric minimum 



scenario the uplift is provided by an anti D3-brane in a throat (in the bulk) of the compactification 
manifold, which gives p = 2 (p = 3). The minimum (j2.7p is then separated from the runaway 
solution by a barrier whose height is set by the gravitino mass. 

Before discussing the dynamics of the system, let us comment on the approximations upon 
which we rely in the remainder of the paper. It is known that (warped) string-inspired construc- 
tions, such as KKLT, receive additional contributions, e.g. a' corrections |21| . open-string loop 
corrections |22| and corrections due to (strong) warping |23y24). Although they may significantly 
alter the form (12. ip of the supergravity action |25| . we do not incorporate these effects in our 
study. Our main concern is to understand whether a gaugino condensate can be integrated in/out 
when the dynamics of a modulus is taken into account. As such, we perform a comparative study 
between the single field behavior and the full system evolution. The corrections mentioned above 
are in this sense somewhat orthogonal to our purpose, and should not affect our conclusions. 
Moreover, let us emphasize again that we focus on the KKLT model as an illustrative example. 
Our observations hold for all moduli stabilization scenarios which rely on gaugino condensation, 
including for instance those without warping. 

2.1 Dynamical stabilization 

The barrier separating the KKLT minimum from the runaway solution at infinity is small com- 
pared to any other natural scales in the problem. This is true, in particular, for the phenomenolog- 
ically favored case of low-energy supersymmetry breaking. At the same time, the scalar potential 
is exponentially steep. One can therefore wonder how likely it is that the modulus will end up in 
the right minimum and not roll off to infinity. This is known as the overshoot problem, and has 
first been recognized in |26| . Dynamically, the higher the energy with which the field starts its 
dynamical evolution, the more friction is needed in order to slow it down so that it finally settles 
in the minimum. 

In |10H14| . it was shown that the presence of a background fluid can help avoiding the over- 
shoot problem by providing the necessary friction. Indeed, there exists a relatively large range of 
initial conditions that lead to successful stabilization. Clearly, this range grows with the initial 
background density, see e.g. Fig. 1 and 2 in [14]. Throughout this paper, we consider the back- 
ground fluid to be composed of radiation. This allows us to work in a model-independent way. If 
the expansion of the universe were instead driven by matter or vacuum energy density, the con- 
densate would acquire a large mass of order with a model dependent coefficient. Moreover, as 
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was shown in |15) . if the background field gives the dominant contribution to SUSY breaking, the 
induced mass term destabihzes the gaugino minimum at a scale H ~ cA, with c <C 1, somewhat 
below the condensation scale (j2.5p . In a radiation dominated universe, the Hubble induced mass 
for the condensate and the volume modulus is small and can be neglected at tree level |27| . For 
instance, a coupling of the form uti*T^', where T is the stress-energy tensor of radiation, vanishes 
at tree level. One might however worry that this is not necessarily the case at one loop. We will 
come back to this point in Section |4] when we consider finite temperature corrections. 

In supergravity, the dynamics of a set of fields with non-canonical kinetic terms is described 
by the equations of motion 

4>' + 3H^' + Ti^,^^^ + K'^djV = 0, (2.9) 

where the Christoffel symbols are T*^ = ' . The Hubble radius is subject to the Friedmann 
constraint: 

= £ki„ + l^ + P7 ' (2-10) 

where C\^\a = Kij(j)^(j>', and the energy density of the background radiation satisfies 

p^ + AHp^ = =^ = p^.inie"''^, (2.11) 

where in the last equality, we have used as time variable the number of e-folds = In i?, with R the 
scale factor (recall that H = R/R). Here and in the following we use a subscript to denote the 
corresponding quantity at the initial time N = 0. Without the background component, any initial 
conditions satisfying T^(c7ini) > ^(<7max) will lead to overshoot of the fields and consequently to a 
phenomenologically unacceptable runaway solution. However, as mentioned before, if the energy 
density in the universe is dominated by radiation, the evolution of the scalar fields is damped, 
and a much larger range of initial conditions leads to modulus stabilization. This effect enters the 
equations of motion (j2.9p via a large Hubble rate H ~ P7/3. 

The explicit equations of motion derived from (I2.9P for the KKLT model are given in Ap- 
pendix [Al Eqs. (lA.SP in the case of the potential ()2.6p can be solved numerically. With this 
purpose, we assume that the field has vanishing initial velocity Ho-, mi = 0. Following |12|I14|. we 
parametrize the initial radiation as 

P^,ini = Y^^V{ar,ini), (2.12) 

where we used the fractional energy density in radiation $7^ = /9^^ini/(3i?^). Note that for 0^ > 1/2 
radiation initially dominates. An example of modulus stabilization is given in Fig. [H where we 
have taken = 0.99. The rest of the parameters is given as in [9l[T2], i.e ^ = 1, a = 0.1 and 
771,3/2 — 1 TeV. The constant Wq and C are deduced from (12. 7p and (12. Sp . 

2.2 Consistency Condition 

As ar rolls down the KKLT potential in the presence of the thermal fluid (j2.1ip . the corresponding 
scale A = Umin = exp [— (aa^ -|- 1) /3] evolves as well. Indeed, only when the modulus is settled 
in its minimum does the gauge coupling of the hidden sector stabilize as well. In order for the 
potential to follow effectively from (12. 6p . we have to ensure that condensation does occur even as 
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Figure 1: Example of initial condition leading to stabilization. The left panel shows the evolution of the 
field with time, in this case cr,., ini — 40. The dashed green line gives the position of the minimum. On the 
right figure is plotted the evolution of A (red, dashed line) and H (blue, solid line) with time. 
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Figure 2: Set of initial conditions {(Tr, ini, cn/mi} leading to stabilization at the minimum of the KKLT 
potential for p-^^ini = 10^^, A = I, a ^ 0.1 and 7713/2 = 1 TeV. The solid line corresponds to the result 
of IH]. The darker region can be excluded since Hiai > Aini. 



the modulus is dynamically evolving. In other words, there is a natural consistency condition: 
the Hubble rate H has to be smaller than the condensation scale A at all times. 

It can be seen from the right plot of Fig. [1] that the scale A for this specific example is indeed 
higher than H and remains as such with time. Actually, because the field starts with vanishing 
velocity, it is clear that if i^ini < Ajni, then the same holds at later times. 

In Fig. [2] we show how the region of successful initial conditions for a found in |12| is modified 
once we require that these points are consistent with having the condensate initially formed. 
Actually, even if the condensate is initially formed, and thus remains so at later times, it is 
legitimate to wonder below which scale it can be safely integrated out. This is why we also show 
the bound Hi^i < 10~^Aini in Fig. [2l Formally, a field can be integrated out once its oscillations 
around its minimum are negligible. Hence, for a large range of initial parameters, consistency 
requires that we integrate the condensate in, and consider the whole system (j2.ip together with 
Wq. In the next section, we let the two fields evolve and study the region of stabilization, putting 
particular emphasis on the way the condensate may be integrated back out. 



6 



3 Dynamics of the gaugino- modulus system 



In this section we discuss the dynamics of the Kahler modulus with the gaugino condensate 
integrated in. For simphcity, in what fohows, we concentrate on the real parts of cr and u, and 
omit subscripts. Dynamics of the imaginary parts shall be treated elsewhere |28| . 

The VY effective potential for gaugino condensate u is given in (j2.ip . the additional uplift 
term in (j2.8p . As the Hubble scale approaches the condensation scale, the evolution of the full 
condensate-modulus system starts to deviate from the KKLT description, which is only a good 
approximation for small u when (12. 3p is valid. This can be seen explicitly from the F-term potential 



-2Wqu^ + (2c + So- + Qclnufu^ + ^(4c + 2a)u^ 



(3.1) 



with X = e~^/^, as before. The quartic term in Vp is proportional to In the limit that 

w <C 1 and Wq small, both the cubic and the term are subdominant, and we recover the VY 
result that Wu = minimizes the potential. The scalar potential is a surface; it exhibits a valley 
of attraction corresponding to the direction nmin(o'). This valley is much steeper in the condensate 
direction than along the modulus direction. Along the u direction there is a barrier separating 
the "good" minimum, from the state at u = 0; in the modulus direction there is the usual barrier 
separating the "good" KKLT minimum from the runaway solution at infinity. 

The equations of motion for a and u momenta in a spatially flat Friedmann-Robertson- Walker 
spacetime are given by ()A.4p . As before, we assume an homogeneous background radiation energy 
density, which evolves in an expanding universe according to (12. lip . 

We integrated numerically the equations of motion for the real fields {o", u} keeping the phases 
fixed at their instantaneous minimum. The free parameters in the model are c, Wq (which de- 
termines the vacuum gravitino mass), and Vt^. In what follows, we choose c = 1/7 and a) 
Wq = 10~^ and b) Wq = 10~^. For these parameters, the minimum with the modulus sta- 
bilized is at fjmin = 1.22(1.68) and uq = 'Umin(o-mm) = 0.01 (0.002) for Wq = 10-^(10"^). 
Demanding a Minkowski vacuum, we obtain the uplifting parameter C ~ 3 • 10"^^ (3 • 10"^'') for 
Wq = 10~^ (10"^). We used relatively large values of Wq, giving rise to high scale supersymmetry 
breaking, because of numerical convenience. Nevertheless, we expect that our qualitative results 
can be straightforwardly extrapolated to lower scales. 

In our study, we take a conservative approach and assume that whenever u jumps the barrier at 
Umax and is driven towards the state at the origin, the gauginos do not condense and, consequently, 
a does not get stabilized. Apart from an understanding of the physics at the minimum « = 0, 
discussing the dynamics of u in this case requires to include quantum effects, such as particle 
production of fields which are light at the special symmetry point at the origin |29) . Particle 
production damp the motion of possibly trapping it in the region \u\ < Umax and leaving our 
results untouched. Another possibility is that the field has enough momentum to go through the 
origin and bounce back over the barrier to end up at itmin- This could open up an additional 
region of preferable initial conditions, which nonetheless will strongly depend on the details of the 
potential near the origin. 



3.1 When can the gaugino condensate safely be integrated out? 

The first issue to address is to determine for which initial energy scales H[^i it is a good approxi- 
mation to integrate out the gaugino condensate and work in the low energy effective KKLT theory, 
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Figure 3: The regions bounded by the curves correspond to initial conditions leading to modulus stabi- 
lization in the KKLT model with c = 1/7 and a) Wq = 10^^ and b) Wq = 10~® as a function of the initial 
abundance Q-y. The thin/red curve is obtained in absence of u and the thick/blue curve is for u integrated 
in and initially at its minimum Umin- Black dotted lines correspond to (i?/A)ini = 1, 10^^, 10^'^. 



and for which initial conditions the dynamics of the full system should be taken into account. To 
answer this question we determined the range of initial conditions for (Jini leading to modulus 
stabilization in the effective KKLT model, with u integrated out. We compared the results with 
the range of fJini leading to stabilization in the full system, where both u and a are kept as dy- 
namical fields, starting with u initially lying in the valley of attraction, i.e. Uini = Umin (fJini). We 
numerically determined this Ui^i, as the analytic approximation (j2.4p breaks down for small Cini 
(which corresponds to large u). In both cases we started with zero momenta for all fields. 

The space of initial conditions {<Tini, ^-y} for which the system evolves to the overall minimum 
is plotted in Fig. [3l for all other initial conditions, u and/or a overshoot. The thin/red curve 
gives the stabilization region for the KKLT potential with u integrated out. The results agree 
with previous findings |12| : provided that the initial energy density is dominated by radiation, 
there is a large range of initial Uini leading to stabilization. Starting with initial values Cini to the 
left of the curve, the modulus picks up too much momentum and overshoots; starting with initial 
values O'ini to the right of the curve, the modulus is too close to the minimum and there is not 
enough time to damp the field, with overshoot as result. Finally, we note that there is another 
small stabilization region for (Jini close to the minimum dtnin (that is not displayed in the plot), 
for which the initial potential energy is less than the height of the barrier Vi^i < l^nax- 

The equivalent curve for the full system is the thick/blue line. As expected, there are large 
discrepancies for small Cini, and thus large Uniin(o"im)i whereas for larger (Jini both systems give the 
same results. In terms of n = {H/A)[^i, the deviations are negligible when this ratio drops below 
n = 10~^, and become large around n = 10"^ — Fig. [3] shows the lines for n = 1, 10~^, 10~^. 
This is almost independent of Wq, i.e. of the vacuum gravitino mass; the reason is that for such 
initial conditions the cubic term in the potential (j3.ip is negligible, and all Wq dependence drops 
out. This allows to restate the results in terms of l^ni as well: for l^ni ^ 10""^ the KKLT and full 
potential behave very similarly, whereas for Vi^i ^ 10"^ they deviate significantly. 

We note that the smaller Wq is, the longer it takes for the fields to roll down to their minima, 
and the longer they are exposed to damping; this explains why less initial radiation is needed 
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for stabilization. Hence, for a TeV gravitino mass, Wq ~ 10"^'' and c = 1/10, it is natural to 
expect the same behavior as in Fig. [3l with an even larger accessible range of $7^. In that case, 
one can further show that stabilization occurs at the unification-compatible point = 1/2. 

As already stressed, the main reason for the breakdown of the KKLT approximation at large 
initial energy densities is due to the way the condensate is integrated out. Formally, the field u 
should be integrated out at the level of the scalar potential. In the KKLT approach, however, 
this is not the case: the field u is sent to zero in the Kahler potential, and replaced by its vev in 
the superpotential. The difference between the two approaches is easy to understand by studying 
the structure of the scalar potential (j2.2p ; it gives rise to extra terms when the first approach 
is followed. The way heavy fields can be integrated out consistently in supergravity has been 
the center of renewed attention in the past months |30H32) . In this paper, we are dealing with 
intermediate energy ranges for which the extra terms in the scalar potential are not negligible. 
They have an impact on the dynamics of the system. Therefore, for small cJini, the solution 
for Umin (|2.4p is no longer a good approximation, and thus higher order terms in u (the u^- 
term in eq. (13. ip ) have to be taken into account. Numerically, we find that the exact and the 
analytical solutions start to deviate more and more as we increase l^ni- Notice that for cJini < 0.2 
(or l^ni ^ 10~^), the nontrivial minimum of VpifT = o"ini, u) along the u direction disappears 
altogether, implying that gaugino condensation cannot occur. This explains the cutoff in the 
stabilization region in both plots. For cJini < 0.2, any Uini leads to overshoot. 

3.2 Initial conditions leading to modulus stabilization 

We now determine the stabilization region in {(Tini? "^ini} space, considering initial conditions with 
both the modulus and the condensate displaced from the instantaneous minimum. Once again 
the fields start their evolution at rest. 

The resuhs for Wq = 10"^ and = 0.95 are shown in Fig.[H The initial conditions {cxini, 6u = 
^ini ifmin(fini)} leading to stabilization lie in the shaded/green area. A large range of modulus 
initial values leads to stabilization 0.2 < cJini ^ 0.6, whereas only a small displacement in 8u is 
allowed. The reason for this is simple: the potential is much steeper in the condensate direction, 
even a comparatively small displacement \6u\ < 2 x 10~^ implies a large increase in Vini, and 
consequently leads to overshoot. 

To address the likelihood of initial conditions leading to stabilization, it is therefore useful to 
consider not the initial field displacement, but rather the corresponding change in energy Vini- 
Given an initial energy l^ni, and random initial conditions for both the modulus and condensate 
(but such that l^ni remains fixed) , how large is the stabilization region in the space of initial con- 
ditions? In Fig. m are plotted the equipotential lines corresponding to Vini = 10~^, 10~^, . . . , 10~®. 
The lines are approximately orthogonal to Su = for small CTini- Considering all initial condi- 
tions along one such line to be equally likely, the stabilization region is very small. Most initial 
conditions lead to overshoot. This is quite a different result than the conclusion one (naively) 
draws from the one-dimensional system with the condensate integrated out. Confining to initial 
conditions to the left of the maximum, a good proportion of initial conditions leads to stabilization. 

The solid/orange line at 5u = in Fig. |4] corresponds to a similar study in the effective KKLT 
model, without the dynamics of u. The points along that line at the left of the shaded/green area 
are initial conditions that appear to lead to stabilization from the one-dimensional standpoint. 
Once the dynamics of the condensate are considered, it is clear that those points are inconsistent. 
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Figure 4: Initial conditions {(Tini,Ui,ii — Wmin(cini)} leading to modulus stabilization for Wq = 10 ^. The 
black dashed lines correspond respectively to Vini = 10^^, 10^^, . . . , 10^^ from left to right. 



4 Finite temperature corrections 

In this section we study the effect of thermal corrections on the evolution of the modulus- 
condensate system. These corrections arise because the gauge kinetic term W"'W"' of the thermal 
sector acquires an expectation value at one loop |33p34j. and because the particles in the thermal 
bath may have a field dependent gauge or Yukawa coupling. As an example, let us consider that 
the background thermal fluid /j^ originates from an SU(iVc) gauge sector with Nf matter multiplets 
(in the fundamental representation) and a modulus-dependent gauge coupling, = \/a. Even 
though the (real) modulus a does not thermalize, its potential receives corrections through the 
gauge coupling of the thermal bath. Hence, at a constant temperature T, the effective potential 
energy of the modulus is 

V{a)+F{g,T), (4.1) 

where V{a) is the zero-temperature potential (j3.ip and the induced free energy is given bjH 

F{g, T) = {ao + a^g^) T\ (4.2) 

where qq = {N^ + 2NcNf - l) and 02 = ^ {N^ - l) (N^ + 3Nf). Note that for a sufficiently 
high temperature Tcrit, the minimum of the potential gets lifted and, if the modulus is already in 
its minimum, it gets destabilized |36y37). 

Similarly, since the trace of the stress-energy tensor no longer vanishes, a mass term for the 
condensate may be generatecQ via a coupling uu* (T^) It is however easy to show that the 
corresponding destabilization temperature is (much) higher than the one affecting the modulus, 
21;rit- We have confirmed numerically that such a mass term can be neglected. 



^Formally, eq. (|4.2|) is valid only in the weak coupling regime. If g ~ additional corrections become 

important. However, explicit computations have shown that the term contains aheady the correct qualitative 
behavior |35| . 

''We thank Oleg Lebedev for drawing our attention to this issue. 
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The existence of a thermal bath affects the moduli dynamics. First, the solution of Friedmann's 
equations for the radiation density becomes 



where r denotes the raticH 02/00 and = Ini? as before. Since enters the Hubble rate H, 
eq. (j4.3p implies that temperature corrections introduce an additional source of friction which 
decreases as the universe expands. As a result, the equations of motion ()2.9p are modified as 
follows 

djV ^ dj{V + F) = djV - I j^djg' , (4.4) 

where we have used the relation between the energy density and the temperature of the thermal 
fluid = —p^ + Tdp^/dT = — 3ao(l + rg^) with the pressure given by p^ = —F{g, T). This 
new contribution to the equations of motion renders manifest the effect of the effective poten- 
tial (j4.ip . Note that irrespective of whether the initial temperature is larger than the temperature 
at which the minimum of the potential gets hfted, Tjni > Tcrit, the modulus can still be stabilized. 
If the friction is large enough, the minimum emerges before the modulus has had time to approach 
it. Some of the initial conditions that lead to stabilization at zero temperature are however erased. 
Thus, the region of admissible conditions reduces under the influence of the thermal bath. 

We note that it is possible that a subset of the thermal degrees of freedom do not have a 
modulus dependent gauge or Yukawa coupling, and only interact gravitationally with the modulus 
sector. If this subset is large, the parameter |r| can be made arbitrarily small, and the results of 
the previous section apply. 



4.1 Initial conditions leading to modulus stabilization 

As before, we consider the KKLT model with the condensate field u integrated in and the pa- 
rameters c = 1/7 and Wq = 10~^. In addition, we take the initial radiation abundance to be 
Q,.y = 0.99 and consider two different configurations with r = — 3/47r^ and r = —0.04 to illustrate 
the dependence on thermal corrections. The first case corresponds to an SU(2) pure Yang-Mills 
thermal sector and is thus the minimal |r| achievable when the contribution (j4.2p is the only 
component in radiation. Fig. [5] shows the set of admissible initial conditions {aini, Su} that lead 
to stabilization. The region surrounded by the dashed line corresponds to the result obtained in 
the previous section for zero temperature (cf. Fig. [4]). Inside this region, we find those admissible 
initial conditions for r = —0.04 (light shaded) and r = — 3/47r^ ~ —0.076 (dark shaded). 

A few comments are in order. The region of stabilization clearly reduces when |r| increases. 
For the SM gauge group SU(3) xSU(2), one finds r ~ —0.1, which leads to a considerable reduction 
of the allowed region. This implies that only for an even more fine-tuned value 0^ 1 does stabi- 
lization occur in a more realistic scenario. Also, the smallest value that the modulus can initially 
take depends on the thermal background. This happens because the higher the potential energy 
of the modulus l^ni is, the faster it rolls down to the position of the minimum. If the temperature 
is not sufficiently low, the modulus finds no minimum and consequently runs away. Similarly, if 

^Our definition of r diflfers by a factor 4tt compared to |14| . This originates from the different dependence of 
the gauge coupling on a. 
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Figure 5: Region of initial conditions leading to a stabilization in the KKLT model. The dashed line 
surrounds all admissible initial values in the absence of thermal corrections. The light-shaded (dark- 
shaded) region corresponds to the allowed parameters in presence of thermal corrections, with r = —0.04 
(r = -3/47r2). The chosen parameters are c = 1/7, Wq = 10"^, Q-y = 0.99. 

a starts too close to the local maximum of the modulus potential, at higher temperatures the 
modulus does not feel the attraction of the minimum and escapes. 

We explore also how the stabilization of the modulus depends on the initial radiation abun- 
dance 0^. In Fig. [6l we present our results including the condensate dynamics (blue/thick curves). 
In the left (right) panel we have set Wq = 10"^ (10~^). The area over the solid (blue/thick) curve is 
the set of stabilization-compatible initial conditions at zero temperature. The area over the dashed 
(blue/thick) curve describes stabilization when thermal effects are included, with r = — 3/47r^. 
As expected, for low initial radiation abundance, the modulus runs away before the temperature 
decreases below Tcrit and, therefore, stabilization is not possible. This sets a bound on the initial 
value ofn^, which in our examples is n^,mm ~ 0.98 (0.92) for Wq = 10"^ (10~^) with r = -S/Att'^. 

We contrast these findings with those results obtained in the traditional way, i.e. in absence of 
the dynamical condensate field u (red/thin curves). Note that the thin and the thick curves only 
coincide in a narrow region close to the stabilization point of a (iTmin ^ 1-22 for Wq = 10~^ and 
fmin ~ 1-68 for Wq = 10~^). This implies that, integrating out the condensate is only valid for 
these points, that is, for initial conditions such that the Hubble rate i^ini is between two and three 
orders of magnitude lower than the condensation scale Aini at the beginning of modulus evolution 
This observation is rather model independent and has to be taken into account in general when 
dealing with moduli stabilization driven by gaugino condensate in an evolving universe. 

5 Conclusions 

We have studied the cosmological evolution of a modulus coupled to a gaugino condensate and its 
consequences for moduli stabilization. Using the KKLT model as an illustrative example, we have 
shown that, in presence of background radiation, most of the region yielding modulus stabilization 
corresponds to points in parameter space for which the Hubble parameter i^ini is close to the scale 
at which the gauginos condense Aini. Once the dynamics of the gaugino pairs is included, we have 
shown that: 

• integrating out the condensate directly in the superpotential and the Kahler potential fails 
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Figure 6: The regions enclosed by the curves correspond to initial conditions leading to modulus stabiliza- 
tion in the KKLT model with c = 1/7 and a) Wq = 10~^ and b) Wo = 10^^ as a function of the initial 
abundance il^. The thin/red curves are obtained in absence of the field u for r = (solid) and r = — 3/47r^ 
(dashed). The thick/blue curves are obtained for u integrated in and initially at its minimum Umin for 
r = (solid) and r = — 3/47r^ (dashed). Black dotted lines correspond to (i7/A)ini = 1, 10~^, 10~^. 



to describe moduli dynamics unless H[^i is about two orders of magnitude smaller than Aini, 

• the dynamics of the condensate sets a stronger bound on the initial conditions leading to 
stabilization. 

The former point implies that if one is to consider that the condensate has been integrated 
out before the modulus starts its evolution, then only the vicinity of the minimum leads to 
stabilization. 

We have also discussed the role of thermal corrections and shown that they do not affect the 
conclusions above. Thermal effects impact on the initial friction needed to damp the modulus 
and stabilize it, constraining further the range of initial conditions that lead to stabilization. 

To conclude, we point out that our observations are model independent and have to be taken 
into account in general when dealing with moduli stabilization driven by gaugino condensate in an 
evolving universe. Also, our results are expected to hold independently of the choice of parameters. 
In particular, they should be valid in settings with a realistic supersymmetry breaking scale. 
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A Equations of motion 

The equations of motion (12. 9p 

4>' + 3H^' + ri,^,j>><j)^ + K'^djV = 0, (A.l) 
can be separated into real and imaginary parts (resp. denoted <j)r and 0j) and read 
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Applied to the KKLT system ()2.6p . we have 



(Tr + 3/i(Tr. (cJr - O-j ) + ^—7^ = , 

'6 oar 

ai + 3Hdi- —drdi + ^^ = . (A.3) 
ar o oai 

In section [3l we aim at comparing the region of stabihzation with the results obtained in |12| . 
However, from then on, we concentrate on the real parts of a and u. Dynamics of the imaginary 
parts shall be treated elsewhere. 

For the exact system (12. ip . the first equation in (|A.2j) becomes 

dr + 'iHdr-- ^ 9 /o i^r^ ~ -^^rdr^ + '"r/3 (20-^9^^^ + Urdur-V) = 0, 

iGr — Uf./6 \ 6 J 

Ur + ?>HUr - i^Ur^ - Urdr) + ^^^^ ~ "^^'/^ (^,9,^^ + 35^,^) = . (A.4) 

lar — uf./6 \ 6 J 

In general, it proves much easier to re-express eqs. (|2.9p in terms of the canonical momenta 
associated to the fields Ilj = 3£kin/9(/'* 

ct>'' = (A.5) 

where prime denotes derivative with respect to (i.e x = Hx'). These equations can then be 
integrated numerically in order to obtain the solutions (ff-^N). 
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